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INTRODUCTION 



This manual describes the Control Data® 3600 FORTRAN library routines. 
Library routines may be called directly by the COMPASS calling sequence 
listed, unless otherwise noted. 

To test the accuracy of the double precision functions, several thousand 
arguments covering the necessary ranges were generated. The functions were 
evaluated, and the answers compared with the corresponding answers ob- 
tained from the CO-OP library routine A3 NBSB DOBLPREC. The latter 
routine does 91-bit interpretive arithmetic. Since, in general, different 
algorithms were used and the answers agreed to an average of 82 bits, the 
accuracy of the double precision functions was considered satisfactory. 

Accuracy of the single precision functions was determined by comparing them 
with the double precision functions in a similar manner. In testing the timing 
of routines, several thousand appropriate arguments were chosen. The 
functions were evaluated for these arguments, and the time measured by the 
TIMEF function. The times quoted within this manual include the jump into 
and out of the routine, but do not include indexing, loads, and stores of the 
main programs. 

Execution times for complex and double precision functions for the Control 
Data 3600 have been computed for several thousand arguments using the 
TIMEF routine. The arguments were chosen to cover uniformly a reasonable 
portion of the allowable ranges. The range from to 7r/2 was timed separately 
for DSIN and DCOS because of a significant difference in execution time. 
In all cases, the time includes the jump into the routine and the jump back, 
but does not include any other instructions in the testing program. 



ABSF 



PURPOSE: To obtain the absolute value of the real or integer argument in the A register. 

COMPASS CALLING SEQUENCE: 
Call by Value (ABSF) 

LDA 
L BRTJ 

L+l SLJ 

00 



Xor I 

($)ABSF, , * or ($)XABSF, ,* 

*+l 

DICT. 



Call by Name (Q8QABSF) 
L BRTJ 

L+l SLJ 

01 
L+2 00 

00 



($)Q8QABSF, , * or ($)Q8QXABSF, ,* 

*+2 

DICT. 

($)X or ($)I 



FORTRAN FUNCTION: ABSF(X) or XABSF(I) 

NORMAL RETURN: Result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 7 locations 

ABSF and XABSF are generally compiled on-line. 
Compiled coding sequence: 

LDA X or I 

+ AJP, PL * + 1 

ROP.XOR A.MZ.A 



ACOSF 

PURPOSE: To compute the single precision arccosine of a single precision argument X. 
COMPASS CALLING SEQUENCE: 



Call by Value (ACOSF) 




LDA 


X 


L BRTJ 


($)ACOSF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QACOSF) 




L BRTJ 


($)Q8QACOSF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: ACOSF(X) 

NORMAL RETURN: Floating point result in the range to n is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN ASIN/ACO 



A = (value) 
ARG GT 1. 



Q = (value) 



If the absolute value of the argument is greater than 1.0, this message is printed 
on the standard output unit, and a normal exit is taken from ACOSF. 



STORAGE: 96 locations 



-34 



ACCURACY: Average 35 bits with a relative error less than 2 maintained throughout the domain 
of the argument. See ASINF. 

TIMING: Average 178 microseconds for 5000 random values of X, -l^X^l, 

MATHEMATICAL METHOD: The identity ACOS(X) = -|- - ASIN(X) is used. 

See ASINF for method description. 



AIMAG 



PURPOSE: To obtain the imaginary part of a complex number Z. 
COMPASS CALLING SEQUENCE: 



BRTJ 


($)AIMAG, , * 


SLJ 


*+2 


01 


DICT. 


00 


($)Z 


00 





FORTRAN FUNCTION: AIMAG(Z) 

NORMAL RETURN: The imaginary part of Z is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 9 locations including REAL and CMPLX 

COMMENT: In a FORTRAN program, Z must be declared type COMPLEX. 

TIMING: 46 microseconds 

AIMAG is generally compiled in-line. 
Compiled coding sequence: 

LDA Z+l 



ALOG10 



L 


BRTJ 


L+l 


SLJ 




01 


L+2 


00 




NOP 



PURPOSE: To compute the common (base 10) logarithm of a floating point argument X. 
COMPASS CALLING SEQUENCE: 

($)ALOG10,,* 

*+2 

DICT. 

($)X 



FORTRAN FUNCTION: ALOG10(X) 

NORMAL RETURN: Returns with the result in the A register. 

ERROR MESSAGE: Since ALOG10 calls LOGF, see LOGF for error messages of arguments less 
than or equal to zero. 

STORAGE: 16 locations 

ACCURACY: 34 bits, relative error increases for X close to 1. 

TIMING: 190 microseconds 

MATHEMATICAL METHOD: 

Log 10 (X) = (Log 10 (e)* Ln (X) 

where Log 10 (e) = Log 10 (2.7182818284) 
= 0.4342944819032 

LOGF is called to compute Ln (X). 



ASINF 



PURPOSE: To compute the single precision arcsine of a single precision argument X. 
COMPASS CALLING SEQUENCE: 



Call by Value (ASINF) 




LDA 


X 


L BRTJ 


($)ASINF, , * 


L+l SLJ 


*+l 


00 


DICT. 


L+2 Return 




Call by Name (Q8QASINF) 




L BRTJ 


($)Q8QASINF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: ASINF(X) 

NORMAL RETURN: Floating point result in the range - — to — is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN ASIN/ACO 
A = (value) Q = (value) 

ARG GT 1. 

If the absolute value of the argument is greater than 1.0, this message is printed 
on the standard output unit and a normal exit is taken from ASINF. 

STORAGE: 96 locations 

ACCURACY: ASINF was compared with a double precision arcsine routine for 10, 000 random num- 
bers between +1 and -1 . The error distribution is as follows : 



No. of bits 
in error 


-4 


-3 


-2 


-1 





1 


2 


3 


4 


Frequency 





24 


293 


2020 


5234 


2072 


339 


15 






TIMING: Average 178 microseconds for 5000 random values of X, -1 s. X« 1. 



ASINF 

MATHEMATICAL METHOD: 

Given an argument | X| between 0. and 1. then 



ASIN (X) = A + BX (k x + X 2 (k 3 + X 2 (k 5 + X 2 (k ? + X 2 (k 9 + X 2 (k n + X 2 (k 13 + X 2 (k 15 X 2 )))))))) 

This is a relative minimax approximation of degree_15. The domain of X(0 to 1) is parti- 
tioned into 4 intervals and the values of A, B, and X are assigned as follows: 



1. =s |xl =£ .5 

2. .5<lxN.866 

3. .866< lxl^.966 

4. .966<|X|<1. 

The octal constants used are: 



X = X, A= 0, B = 1. 

X = 2X 2 -1, A = -J-, B=.5 

X = 8X 4 -8X 2 + 1, A = ^ L , B = .25 



r i-txn i/2 , 7T 

L 2 J ,A " 2 ■ 



ki 


= 2000 


7777 


7777 


7777 


k 3 


= 1775 


5252 


5253 


1306 


k 5 


= 1774 


4631 


4506 


5155 


k 7 


= 1773 


5556 


7611 


7211 


k 9 


= 1772 


7550 


1757 


0610 


kn 


= 1772 


6246 


7740 


7323 


kl3 


= 1770 


7313 


2650 


2552 


kl5 


= 1773 


4362 


4544 


0473 



ATANF 



PURPOSE: To compute the floating point arctangent of the floating-point argument X. 
COMPASS CALLING SEQUENCES: 



Call by Value (ATANF): 
LDA 

L BRTJ 

L+l SLJ 

00 



X 

($) ATANF, , * 

*+l 

DICT. 



Call by Name (Q 8QATANF) : 



L 


BRTJ 


($)Q8QATANF, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 
NOP 


($)X 



FORTRAN FUNCTION: ATANF(X) 

NORMAL RETURN: Floating point result in the range - — to — is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 64 locations 

ACCURACY: 34 bits 

TIMING: 125 microseconds 

MATHEMATICAL METHOD: The Maehly formula for arctangents is used for most arguments. 

(Control Data Publication No. 516) 

3 



E x t 



tan (X) = C + t D. + t D. + -5- 

n 1 t 2 + D + 







t 2 + D 3 



.0195 < |X| < V2 - 1 

V2 - 1 < |X| < V2 H 1 

V2 + 1 < |X| < 3076 



t = X, C = 
IXl-1 



t = 



t = 



c = — 
|X| 4t ' C 4 



IXI 



«=--§■ 



The formula is not used outside the above ranges. Instead, the following equations are used: 

=s |x| s .0195 
3076 s |X| ^ =° tan X = ± ( - 



1 X 3 x 5 

tan" 1 X = X -^- + ^- 



,2 1X1/ 



ATAN2 



PURPOSE: To compute the angle when the tangent is given in terms of coordinates, X and Y. 



COMPASS CALLING SEQUENCE: 
L BRTJ 

L+l SLJ 

02 
L+2 00 

00 



($)ATAN2, , * 

*+2 

DICT. 

($)Y 

($)X 



FORTRAN FUNCTION: ATAN2 (Y,X) 

NORMAL RETURN: The floating point result in the range -ir to ir of the angle in radians is left in the 
A register. 

ERROR MESSAGE: ATAN (0/0) Undefined. 

STORAGE: 44 locations including CANG and CANGQ8Q 

ACCURACY: 33 bits 

TIMING: 212 microseconds average 

MATHEMATICAL METHOD: ATAN2 establishes the correct quadrant, and calls ATANF. 



CABS 



PURPOSE: To compute the magnitude or modulus of a complex number; real part in Z, imaginary 
part in Z + 1. 

COMPASS CALLING SEQUENCE: 

Call by Value (CMAGQ8Q) (callable from COMPASS only) 

DLDA Z 

L BRTJ ($)CMAGQ8Q, , * 

L+l SLJ *+l 

00 DICT. 



Call by Name (CABS) 




L BRTJ 


($)CABS, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)Z 


NOP 






FORTRAN FUNCTION: CABS (Z) 

NORMAL RETURN: The magnitude is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 18 locations 

ACCURACY: 34-35 bits 

TIMING: 203 microseconds 

MATHEMATICAL METHOD: For the complex number Z, where Z = X + iY 

CABS (Z) = | Z | = SQRTF (X*X+Y*Y) 
CABS (Z) is evaluated as: 
|x| *SQRTF (1. + (Y/X)**2) if |x| > | y| 
or as |Y| *SQRTF (1. + (X/Y)**2) if |y| > |x| 

NOTE: Exponent overflow may occur if |x| or | Y| >6.35*10 307 . If 

overflow occurs, no error indication is given, and CABS (Z) will be 
set to a very small number (of the order of 10~ 308 ). 



CANG 



PURPOSE: To compute the argument or angle of a complex number, such that - it < 6 < ir, 9= argument 
of Z. The real part is in Z, the imaginary part in Z + 1. 

COMPASS CALLING SEQUENCE: 

Call by Value (CANGQ8Q) (Callable from COMPASS only) 

DLDA Z 

L BRTJ ($)CANGQ8Q, , * 

L+l SLJ *+l 

00 DICT. 



Call by Name (CANG) 




L BRTJ 


($)CANG, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)Z 


NOP 





FORTRAN FUNCTION: CANG (Z) 

NORMAL RETURN: The angle in radians is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 44 locations including ATAN2 TIMING: 213 microseconds 

ACCURACY: Accuracy depends on ATANF and the arguments used. 



MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY = R*e ie X, Y, R are real 

Argument of Z = B = Tan _1 (Y/X) such that - it < 6 < n 

1.1 X>o, = Tan- 1 (Y/X) 

1.2 X< 0, Y = 0, e = it 

1.3 X < 0, Y > 0, 6 = Tan -1 (Y/X) + tt 

1.4 X < 0, Y < e = Tan -1 (Y/X) - n 

1.5 X= 

1.51 Y = 6 is undefined, set 0-0, exit 

1.52 Y > 6 = it/2 

1.53 Y < 6 = -tt/2 

where ATANF is used to compute Tan -1 (Y/X). 



JO 



CATAN 



PURPOSE: To compute the arctangent of a complex number; real part in Z, imaginary part in Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($) CATAN , , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CATAN (Z) 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See LOGF. 

STORAGE: 33 locations TIMING: 606 microseconds 

ACCURACY: Accuracy depends on LOGF, ATANF, and the arguments used. 

MATHEMATICAL METHOD: 

For the complex number" Z 

u = real part of ( — J u and w are real. 

♦ * / i + Z \ / i + Z \ 

w = imaginary part of I . I : i.e., 1 , 1 = u + iw 

Tan~ 1 (Z) = \ Ln {jZ^j = ^ Ln (u + iw) = | (Ln(R*e^) ) 



T«n-l(Z)-|m(R) + i*i*--f + i!f® 

_i w 
6 = Tan (— )such that - w< 6 < ir 
u 

R = (u 2 +W 2 ) 1 / 2 

LOGF and CANG are called to perform part of this computation. 



?! 



ccos 



PURPOSE: To compute the cosine of a complex number; real part in Z, imaginary part of Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)CCOS, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CCOS (Z) 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See SINF, COSF, and/or EXPF. 

STORAGE: 36 locations including CSIN TIMING: 624 microseconds 

ACCURACY: Depends on SINF, COSF, and EXPF, and the arguments used. 

MATHEMATICAL METHOD: For the complex number Z 

COS (Z) = SIN (it/2 - Z) 
See mathematical method of CSIN to compute SIN (ir/2 - Z). 



12 



CEXP 



PURPOSE: To compute the natural exponential of a complex number; the real part in Z, the imaginary 
part in Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)CEXP, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CEXP (Z) 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See EXPF, SINF and COSF. 

STORAGE: 21 locations TIMING: 489 microseconds 

ACCURACY: Depends on EXPF, SINF, COSF and the arguments used. 

MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY, X, Y are real 

e Z = e X+iY = e X* e iY = e X (COS(Y) + i SIN (Y) ) 

e Z = e X COS(Y) + ie X SIN(Y) 

EXPF, COSF, and SINF are called to perform this computation. 
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CLOG 



PURPOSE: To compute the natural logarithm of a complex number; the real part in Z, the imaginary 
part in Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)CLOG, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CLOG (Z) 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See LOGF. 

STORAGE: 20 locations 

ACCURACY: Depends on LOGF, SQRTF, ATANF, and the arguments used. 

TIMING: 550 microseconds 

MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY = R*e ie X, Y and R are real 

R = (1 + (Y/X) 2 )l/2 

6 = Tan -1 (Y/X), such that - ir < < n 

Ln(Z) = Ln(R*e ie ) = Ln(R) + 10 

CABS, LOGF, and CANG are called to compute R, 
Ln(R), and 8 , respectively. 



U 



CMPLX 



PURPOSE: To create a complex number from two real numbers, X and Y. 



COMPASS CALLING SEQUENCE: 
L BRTJ 

L+l SLJ 

02 

L+2 00 

00 



($) CMPLX, , * 

*+2 
DICT. 
($)X 
($)Y 



FORTRAN FUNCTION: CMPLX (X,Y) 

NORMAL RETURN: X is left in the A register, Y is left in the Q register. 

ERROR MESSAGE: None 

STORAGE: 9 locations including REAL and AIMAG 

TIMING: 54 microseconds 

CMPLX is generally compiled in-line. 
Compiled coding sequence: 

LDA X 

LDQ Y 
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CONJG 



PURPOSE: To compute the conjugate of a complex number Z; the real part is in Z, the imaginary 
part in Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($) CONJG, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CONJG(Z) 

NORMAL RETURN: The real part of the conjugate is left in the A register; the imaginary part in the 
Q register. 

ERROR MESSAGE: None 

STORAGE: 4 locations TIMING: 50 microseconds 

ACCURACY: Exact 

MATHEMATICAL METHOD: For the complex number Z 

Z=X+i*Y X and Y are real 

Conjugate ofZ=X-i*Y 

CONJG is generally compiled in-line. 
Compiled coding sequence: 

LDA Z 

LQC Z+l 
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COSF 



PURPOSE: To compute the cosine of a floating point argument X in radians. 
COMPASS CALLING SEQUENCE: 



Call by Value (COSF) 




LDA 


X 


L BRTJ 


($)COSF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QCOSF) 




L BRTJ 


($)Q8QCOSF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: COSF(X) 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN SIN/COS 

A = (value) Q = (value) 

ARG GT 2**36. 

If the absolute value of the argument exceeds 2 36 , this message is printed on the 
standard output unit, and the absolute value of the argument is returned in the A 
register. 

STORAGE: 13=9 locations 

ACCURACY: Average accuracy in the range (- -|- , -|-) exceeds 35 bits; the worst known case is 34 bits. 

When |x| > — , the relative error is less than 2 u . 

TIMING: 122 microseconds 

MATHEMATICAL METHOD: 

The method of partitioned polynomials, described in Control Data Technical Report 52 is used. 
Briefly, for |x| > approximately 40°, a polynomial in x is evaluated; for approximately 40° 
Z_|x |Z.90°, a polynomial in (-j- - x) is evaluated; in both cases, the actual angle determines 
how many terms in each polynomial are evaluated. 
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COTF 



PURPOSE: To compute the cotangent of the floating point argument X in radius. 

COMPASS CALLING SEQUENCE: 
Call by value (COTF) 

LDA X 

L BRTJ ($)COTF, , * 

L+l SLJ *+l 

00 DICT. 



Call by name (Q8QCOTF) 




L BRTJ 


($)Q8QCOTF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: COTF (X) 

NORMAL RETURN: Floating point result is in the A register 

ERROR MESSAGE: IF the argument is zero: 

ERROR DETECTED IN COTF 

A = (value) Q = (value) 

ARG = 

STORAGE: 19 locations 

ACCURACY: See TANF 

TIMING: 20 microseconds + TANF timing 

MATHEMATICAL METHOD: The cotangent is evaluated by using the equation 

COTF(X) = TANF(tt/2-X) 
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CSIN 



PURPOSE: To compute the sine of a complex number; the real part is in Z, the imaginary part in 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)CSIN, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CSIN (Z) 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q registe 

ERROR MESSAGE: See SINF, COSF and/or EXPF. 

STORAGE: 36 locations including CCOS. 

ACCURACY: Accuracy depends on SINF, COSF, and/or EXPF 

TIMING: 550 microseconds 

MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY X, Y are real 

Y + -Y ./ Y -Y \ 
CSIN (Z) = 6 - 6 * SIN PC) + ^-y 9 j* COS (X) 
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CSQRT 



PURPOSE: To compute the square root of a complex number; the real part is in Z, the imaginary part 
in Z + 1. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($) CSQRT, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)Z 

NOP 



FORTRAN FUNCTION: CSQRT (Z) 

NORMAL RETURN: The real part is less than or equal to zero. It is left in the A register and the 
imaginary part in the Q register. 

ERROR MESSAGE: None 

STORAGE: 20 locations 

ACCURACY: Approximately 35 bits 

TIMING: 348 microseconds 

MATHEMATICAL METHOD: For the complex number X, let Z = X + iY, 

W = \f~Z~ = U + iV 



-Ic 



Compute C = \j C |x| + \Jx 2 + Y 2 /2 

Then U = CandV=DifX>0 

U = |D| and V = sign (Y)*C if X < 

NOTE: Exponent overflow may occur if |x| or |y| >3.73*10 307 , 
CSQRT (Z) will be set to a number of the order (10 154 , 10 154 ). 
If overflow occurs and if | Z| < 8.98 * lO 30 ?, CSQRT (Z) is 
set to (0,0). 
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CUBERTF 



PURPOSE: To compute the cube root of a real argument X. 
COMPASS CALLING SEQUENCE: 
Call by Value (CUBERTF) 

LDA X 



L BRTJ 


($)CUBERTF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QCUBER) 




L BRTJ 


($)Q8QCUBER, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: CUBERTF (X) 

NORMAL RETURN: The floating point result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 61 locations 

ACCURACY: 35 bits average. Most cases are exact; worst known error is 2 bits 



TIMING: 196 microseconds 



MATHEMATICAL METHOD: 



x l/3 = (2 3n+k * ^1/3 where l/2 < f < lf n and k are like-signed 
integers, and k = 0, ± 1, ± 2 . An initial approximation to f 1 ' 3 is 
given f = A + Bf + Cf 2 , where A = .6937488349, B = .4944862032, 
and C = -.1886870774. Two Newtonian iterations of the form 
f i+1 (2/3) (fj + f/2^ 2 ) are done, the first in floating point, and the 
second using fractional arithmetic to avoid a rounding error. f 2 is 
fractionally multiplied by 2 k , the new exponent is packed, and the 
result returned in the A register. 
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DABS 



PURPOSE: To obtain the absolute value of the double precision argument X. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)DABS, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DABS (X) 

NORMAL RETURN: Result is left in the AQ register. 

ERROR ME SSAGE : None 

STORAGE: 5 location 

TIMING: 40 microseconds 

DABS is generally compiled in-line. 
Compiled coding sequence: 

LDA X 

AJP.PL * + 2 

+ ROP.XOR A.MZ.A 

ROP.XOR Q.MZ.Q 
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DATAN 



PURPOSE: To compute the double precision arctangent of a double precision argument X. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($) DATAN, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DATAN (X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGES: None 

STORAGE: 107 locations TIMING: 898 microseconds 

ACCURACY: Average 83 bits, maintains at least 82 bits accuracy. 

MATHEMATICAL METHOD: 1 

For < x < 0.057, a continued fraction expansion is used: 

x/ M x 2 / 

1) arctan x = r^ + 2 y—' 

/D o i=1 D i 

where M=6 for < |x| < 0.024 and M=7 for 0.024 < |x| < 0.057 

For 0.057 < | x | , the identity used is: 

y + Yo 

2) arctan y = - arctan y Q + arctan — 

7T 

For 0.057 < |x| < l,lety=x. For 1 < | x |, let y=l/x, then arctan x = -- arctan y. 

If < | y | < -fi -1, let y = -tan (tt/16). 
If V¥ -1 < | y | < 1, let y = -tan (3tt/16). 
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DATAN 



y + y 

Let W = ~ . Then 

i-yy 

w/ 13 w 2 / 

arctan W = y r- / + 2 /—> 
/bo 1=1 ZbT 



w 



ta o + 



W^ 



b 1 + 



w 



,2 



b 2 + 



W 2 



^ee "Some Basic 1604 Mathematical Sub-Routines, with Analysis of Chebyshev Polynomials Modified 
Taylor Series, and Continued Fractions", Control Data Corporation publication 61 for a discussion of 
the subdivisions of the range; and Ralston and Wilf's "Mathematical Method For Digital Computers" 
p. 30, for a discussion of continued fractions; also D. Teichroew's "Use of Continued Fractions in 
High Speed Computers", NBS report 1514, NBS Project 3011-60-0002, March 1952. 



The octal constants are: 



b Q = 2001400000000000 



h 1 = 2002600000000000 



b 2 = 2001500000000000 



bg = 2002616161616161 



b 4 = 2001504000000000 



b 5 = 2002620376671352 



b 6 = 2001505000000000 



b ? = 2002621173230435 



b 8 = 2001505310000213 



2002621444355714 



J 10 



2001505457754524 



b n = 2002622437536621 



b 12 = 2001533551505605 



0000000000000000 
0000000000000000 
0000000000000002 
6161616161627261 
0000000015365211 
3035475047030314 
0013510547076322 
1112472750662701 
2616107512023450 
6265443712724534 
5260666322154061 
6257645425724336 
0740405770000000 
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DATAN2 



PURPOSE: To compute the angle when the tangent is given in terms of the double precision argu- 
ments, X and Y. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)DATAN2, , * 

L+l SLJ *+2 

02 DICT. 

L+2 00 ($)Y 

00 ($)X 

FORTRAN FUNCTION: DATAN2 (Y, X) 

NORMAL RETURN: The double precision result of the angle in radians is left in the AQ register. 

ERROR MESSAGE: ARCTAN(0/0) UNDEFINED. 

Written on the standard output unit; a normal exit is taken from DATAN2 when 
both arguments are zero since the arctangent (0/0) is undefined. 

STORAGE: 57 locations 

ACCURACY: 82 bits 

MATHEMATICAL METHOD: To compute 

TAN -1 (Y/X) such that -it < 8 S n 

1.1 X > 0, = Tan -1 (Y/X) 

1.2 X < 0, Y= 0, 8 = w 

1.3 X < 0, Y > 0, 8 = Tan _1 (Y/X) + ■k 
1.4X < 0, Y < 0, 8 = Tan- 1 (Y/X) - it 
1.5 X= 

1.51 Y = 0, 6 is undefined, take error exit 

1.52 Y > 0, 8 = 7r/2 

1.53 Y < 0, 8 = tt/2 

DA TAN is used to compute Tan -1 (Y/X). 
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DBLE 



PURPOSE : To convert a single precision floating point argument X, to double precision. 



COMPASS CALLING SEQUENCE: 




Call by 


value 


(DBLE) 




L 




LDA 


X 


L+l 




BRTJ 


($)DBLE, , * 


L+2 




SLJ 


*+l 


Call by 


name 


(Q8QDBLE) 




L 




BRTJ 


<$)Q8QDBLE„* 


L+l 




SLJ 


*+2 






01 


DICT. 


L+2 




00 


($)X 






NOP 






FORTRAN FUNCTION: DBLE (X) 

NORMAL RETURN: The double precision result is left in the AQ register. 

ERROR MESSAGE: None 

STORAGE: 24 locations including SNGL 

ACCURACY: 84 bits 

DBLE is generally compiled in-line. 
Compiled coding sequence: 

LDA X 

LRS 48 

LLS 48 
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DCOS 



PURPOSE: To compute the double precision cosine of a double precision argument X in radians. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)DCOS, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DCOS (X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN DSIN/DCO 

A = (value) Q = (value) 

X GR THAN 2**83. 

When the argument is greater than 2 83 , this message is written on standard 
output and a normal exit is taken from DCOS. 

STORAGE: 122 locations including DSIN 

ACCURACY: Average 83 bits; worst known relative error is 4 x 2~ 84 . 

TIMING: DCOS (0 to tt/2) 518 microseconds 

DCOS (ir/2 and up) 645 microseconds 

MATHEMATICAL METHOD: The cosine is computed using the identity Cos R=Sin(f - R). 

Refer to method in the DSIN description. 
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DCUBRT 

PURPOSE: To compute the double precision cube root of a double precision argument X. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($) DCUBRT, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DCUBRT (X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGE: None 

STORAGE: 55 locations 

ACCURACY: Average 83 bits 

TIMING: 331 Microseconds 



_ o3n 



MATHEMATICAL METHOD: For |X| = 2 



1/8 < f < 1 



The coefficients for the united Chebyshev approximation are those 
generated by Stoer's version of the Remez algorithm 1 : 



Yo = A+B*f + — t 



A 
B 
C 
D 



0.8898 7969 

0.2438 0614 

0.1698 8794 

0.2795 9498 



The Newtonian iteration is of the form 



Y i+1 _ 2 



3f 



f 



2 
4Y, + „ 
1 Yi/2 



Yo and Yi are evaluated with single precision; Y 2 with double precision 

Then 3 Vx = 2 n * Y 2 * sign (X) 

NOTE: Refer to CO-OP routine B4 ANL ANL B453 DCUBRT. 



1 J. Stoer, "A Direct Method for Chebyshev Approximation by Rational Functions, JACM , Vol. 11, No. 1 
(Jan. 1964), pp. 59-69. 
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DEXP 



PURPOSE: To compute the double precision natural exponential of a double precision argument X. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)DEXP, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DEXP (X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN DEXP 

A = (value) Q = (value) 

ARG GR THAN 709. 

When argument is greater than 709.089, EXP(X) will exceed the largest floating 
point number. This message is printed on standard output and a normal exit is 
taken from DEXP. 

STORAGE: 88 locations TIMING: 482 microseconds 

ACCURACY: Average 82 bits; worst known relative error is 4 x 2 -84 . 

MATHEMATICAL METHOD: X = M Ln2 + R M = nearest integer to X/Ln2. 

|R| < Ln2/2 

X --» R 

= 2M + 



Ln2/2 Ln2/2 

e x = 2 M e R 



v. R p+s 

where: e = 

ir — o 



A Q + R 2 (A 2 + R 2 (A 4 + R 2 (A 6 + R 2 * A g ) ) ) 
R(A x + R 2 (A g + R 2 (A 5 + R 2 * A ? ) ) ) 
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DEXP 



This rational approximation is a Pade approximation, see Ralston and 
Wilf, "Mathematical Methods for Digital Computers", p. 13. 

This approximation is also a special case of the Modified Taylor Series 
of Obrechkoff, see "Some Basic 1604 Mathematical Subroutines, With 
Analysis of Chebyshev Polynomials, Modified Taylor Series, and 
Continued Fractions", Control Data pub. 61, p. AA19. 



The octal constants used are: 

A Q = 20357567 02100000 



A x = 20345270 14252363 



A 2 = 20316736 51040243 



A„ = 20265374 55112340 



A. = 20225667 32357466 
4 

A c = 20164251 50626015 
o 



A e = 60114273 66467055 
o 



A ? = 20035421 



A = 17736644 



67776704 
04650356 



00000000 
22265012 
74724732 
05736560 
45736172 
02447135 
13742110 
21366464 
37521044 



00000000 
00634545 
70047333 
14746154 
43317057 
14134622 
25554071 
13166547 
40013400 
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DIMF 



PURPOSE: To determine the difference between two floating point numbers, X* and X„. 
COMPASS CALLING SEQUENCE: 



Call by Value (DIMF) 




LDA 


X l 


LDQ 


X 2 


L BRTJ 


($)DIMF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QDIMF) 




L BRTJ 


($)Q8QDIMF, , * 


L+l SLJ 


*+2 


00 


DICT. 


L+2 00 


($)x x 


00 


($)x 2 



FORTRAN FUNCTION: DIMF(X ,X ) 

X it 

NORMAL RETURN: If Xj > Xg, the floating point difference is left in the A register. 
If X, — X„, zero is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN DIMF 

A = (value) Q = (value) 

OVERFLOW ERROR. 

This message is written on the standard output unit and a normal exit is taken 
when the exponent of the difference, X.. - X 2 , exceeds 308. 

STORAGE: 33 locations 

ACCURACY: Depends on relative magnitude of X. and X 2 

DIMF is generally compiled in-line. 
Compiled coding sequence: 

LDA XI 

FSB X2 

+ AJP, PL *+l 

ENA 



DLOG 



PURPOSE: To compute the double precision natural logarithm of a double precision floating 
point argument X. 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)DLOG, ,* 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DLOG (X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN DLOG 

A = (value) Q = (value) 

ARG = 0/NEG. 

This message is written on standard output and a normal exit is taken from 
DLOG when the argument is zero or negative. 

STORAGE: 124 locations TIMING: 442 microseconds 

-82 
ACCURACY: Average 82 bits; close to 1.0 (from .88 to 1.2), an absolute error of 2 °* is maintained 

but the relative error may be larger in this range. 
MATHEMATICAL METHOD: X = 2* 1 * f 1/2 < f < 1 

Ln f = (A + At 3 + At 5 + At 7 + A / + A /) - a * Ln2 

t - (f - B) / (f + B) 
Ln X = Ln f + n*Ln2 



The range is subdivided into eight intervals yielding: 

m (I-l)/8 /ixC/ 8 * 2*1-1 /1\ 

(I) ^(2) -thenar— and B=(-) 



for 



This makes use of logarithmic subdivision with n=8, see Maehly, Hans J. 
"Approximations for Control Data 1604". The approximation above is a 
Taylor power series of degree 15 telescoped into an 11th degree series 
by the use of Chebyshev polynomials. See "Some Basic 1604 Mathemat- 
ical Subroutines with analysis of Chebyshev Polynomials, Modified Taylor 
Series, and Continued Fractions", Control Data pub. 61, p. AA20. Also 
see Hammings "Numerical Analysis". 
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DLOG 



The octal constants used are: 

A x = 20017777 77777777 77777777 77777703 

A 2 = 20005252 52525252 52525253 05377105 

A g = 17766314 63146314 63113252 20625157 

A 4 = 17764444 44444452 33021367 57217246 

A 5 = 17757070 70570326 44106520 16001122 

A 6 = 17755646 70115712 66255342 54525145 
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DLOG10 



PURPOSE: To compute the common (base 10) logarithm of a double precision argument X. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)DLOG10, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: DLOG10 (X) 

NORMAL RETURN: The double precision result is left in the AQ register. 

ERROR MESSAGE: Since DLOG10 calls DLOG, see DLOG for error messages for arguments less 
than or equal to zero. 

STORAGE: 17 locations 

82 
ACCURACY: Average 82 bits, close to 1. (. 88 to 1. 2) an absolute error of 2 is maintained but 

the relative error may be larger in this range. 

MATHEMATICAL METHOD: LOG 10 (X) = (Log 10 (e) ) * Log e (X) 

Log 10 (e) = Logio (2.71828...) 

= 0.4342944819032518276511289. . . 
DLOG is called to compute Log e (X). 
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DMAX1 



PURPOSE: To determine maximum value of the double precision arguments Xj, X2 X n . 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)DMAX1,,* 

L+l SLJ *+m 

n DICT. 

L+2 00 ($)X 1 

00 ($)X 2 



00 ($)X n 

m = (n + l)/2 + 1, n = number of parameters 
FORTRAN FUNCTION: DMAX1 (X l5 X 2 ,..., X„) 

NORMAL RETURN: The largest double precision argument is left in the AQ register. 
ERROR MESSAGE: None 
STORAGE: 47 locations for DMAX1 and DMIN1 
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DMIN1 



PURPOSE: To determine the minimum value of the double precision arguments X 1 , X 2 , . . . , XqI 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)DMIN1,,* 

L+l SLJ *+m 

n DICT. 

L+2 00 ($)X X 

00 ($)X 2 



00 (ftXn 

m + (n + l)/2 + 1 , n = number of parameters 

FORTRAN FUNCTION: DMIN1 (X x , X 2 X^ 

NORMAL RETURN: Returns with the minimum double precision argument in the AQ register. 

ERROR MESSAGE: None 

STORAGE: 47 locations for DMIN1 and DMAX1 
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DMOD 



PURPOSE: To compute the value of double precision argument X^ modulo X2 


COMPASS CALLING SEQUENCE: 




L 


BRTJ 


($)DMOD, , * 


L+l 


SLJ 


*+2 




02 


DICT. 


L+2 


00 


($)X X 




00 


($)X 2 



FORTRAN FUNCTION: DMOD (X 1> X 2 ) 

NORMAL RETURN: The value of X± modulo X 2 is left in the AQ register. The result is negative 
only if Xj and X 2 have unlike signs. 

ERROR MESSAGE: ARG2 = 0, ARGl/0 UNDEFINED. 

Written on the standard output unit if Xg = and a normal exit is taken. 

STORAGE: 55 locations 

ACCURACY: Depends on the size of the arguments and the subtraction, Xj - Tx^/^l * X 2 . 

MATHEMATICAL METHOD: The result is computed as X x - [Xj/Xg] * X 2 , where only the integer 

part of the bracketed quantity is used. 
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DSIGN 



PURPOSE : To transfer the sign of double precision argument X 2 to the absolute value of double 
precision argument Xj. 



COMPASS CALLING SEQUENCE: 




L 


BRTJ 


($) DSIGN, , * 


L+l 


SLJ 


*+2 




02 


DICT. 


L+2 


00 


($)X X 




00 


($)X 2 



FORTRAN FUNCTION: DSIGNpC^Xg) 

NORMAL RETURN: Result is returned in the AQ register. 

ERROR MESSAGE: None 

STORAGE: 19 locations 

ACCURACY: Exact 
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DSIN 



PURPOSE: To compute the double precision sine of a double precision argument X in radians. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)DSIN, ,* 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)X 

NOP 



FORTRAN FUNCTION: DSIN(X) 

NORMAL RETURN: The result is left in A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN DSIN/DCO 
A = (value) Q = (value) 

X GR THAN 2**83. 

83 
When the argument is greater than 2 , this message is written on standard 

output and a normal exit is taken from DSIN. 

STORAGE: 122 locations including DC OS 

-84 
ACCURACY: Average 83 bits, worst known relative error is 4 * 2 

TIMING: DSIN (0 to ?r/2) 529 microseconds; DSIN (7r/2 and up) 652 microseconds 

MATHEMATICAL METHOD: The sine is computed by using the approximation: 

SIN (R) = R(S + R 2 (S + R 2 (S + R 2 (S + R 2 (S + R 2 (S + R 2 

_L £i O t: O O 

(S ? + R 2 (S g + R 2 (S 9 + R 2 (S 1Q + R 2 (S n + R 2 * S^) )))))))))) 

where < | R| < n/2. 

This is a power series of order 27 telescoped into a series of order 23. 
The telescoping is done through the use of Chebyshev Polynomials. See 
"Some Basic 1604 Mathematical Subroutines with analysis of Chebyshev 
Polynomials, Modified Taylor Series, and Continued Fractions", Control 
Data pub. 61, p. AA20. Also see Hammings "Numerical Analysis". 
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The octal constants used are: 



20014000 00000000 



S = 60022525 25252525 



S = 17714210 



42104210 
60141377 13771377 



S c = 17555616 74351253 
5 

S = 60311214 67246005 
o 

S ? = 17375411 10604724 

S„ = 60501214 01406156 



S n = 17176251 30731453 

S = 60703205 55051622 

S = 16765614 61541561 

S = 61122233 10216717 



00000000 
25252525 
42104210 
13771377 
30714671 
27316212 
13221602 
31461116 
04614442 
44471070 
06010527 
44131274 



00000000 
25252544 
42076135 
20310533 
43064642 
12561732 
12742761 
43664522 
22665460 
34065063 
52104350 
67045655 
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DSQRT 



PURPOSE: To compute the double precision square root of a double precision argument X. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($) DSQRT, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)X 

NOP 



FORTRAN FUNCTION: DSQRT(X) 

NORMAL RETURN: The result is left in the A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN DSQRT 

A = (value) Q = (value) 

NEG ARG. 

When the argument is negative this message is written on standard output and a 
normal exit is taken from DSQRT. 

STORAGE: 37 locations TIMING: 212 microseconds 

ACCURACY: Average 83 bits; worst known relative error is 2 x 2 84 

MATHEMATICAL METHOD: X = 2 2N X, - < X < 2 

Li 

Vx = 2 N * VX 

* X is computed using the Newtonian approximation. 

The initial value is a Chebyshev approximation due to R. F. King 1 : 

Y Q = ^ + 0.464841464 

Two single precision iterations follow: 

Yi = V2 (y { _ 1+ |_ J i = l,2 

Two double precision iterations follow, using the same expression 
with i = 3,4, and giving 

Y 4 ~ VI 
NOTE: Refer to CO-OP routine B4 ANL ANL B451 DSQRT. 



1 Library Subroutine ANL B451 DSQRT Argonnc Natl. Laboratory 
Argonne, Illinois 
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EXPF 



PURPOSE: To compute the floating point exponential of a floating point argument X. 
COMPASS CALLING SEQUENCE: 



Call by Value (EXPF) 




LDA 


X 


L BRTJ 


($)EXPF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QEXPF) 




L BRTJ 


($)Q8QEXPF, ,* 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 





FORTRAN FUNCTION: EXPF(X) 

NORMAL RETURN: The floating point result is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN EXPF 
A = (value) Q = (value) 
ARG GT 709. 

If the argument is greater than 709.089 this message is written on the standard 
output unit and a normal exit is taken from EXPF. 

STORAGE: 42 locations 

ACCURACY: 34 bits, accuracy decreases as arguments get large 

TIMING: 140 microseconds 
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MATHEMATICAL METHOD: 



X = M* Ln2+R*Ln2 M = nearest integer to X/Ln2 

I'RI s 1/2 

X M . RLn2 
e = 2 * e 

OT „ u P + u P/Ln2+R 
Let u = RLn2 ; e = ■ 



P - u P/Ln2-R 

P = A+u 2 (B + r 2") 

V D+u z / 

P/Ln2=A/Ln2+R 2 */B*Ln2+^/ Ln2 ^ 



D/Ln2 2 +R 2 



In octal: 

A/Ln2 = 2002561250731226 
B*Ln2 = 1773433546344025 
C/Ln2 = 2003704560677012 
D/(Ln2*Ln2) = 2007535620560567 
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FLOATF 

PURPOSE: To convert a fixed point number I to floating point. 
COMPASS CALLING SEQUENCE: 



Call by Value (FLOATF) 




LDA 


I 


L BRTJ 


($)FLOATF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QFLOAT) 




L BRTJ 


($)Q8QFLOAT , , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)I 


NOP 






FORTRAN FUNCTION: FLOAT F(I) 

NORMAL RETURN: The floating point result is left in the A register. 
ERROR MESSAGE: None 
STORAGE: 12 locations 

FLOATF is generally compiled in-line. 
Compiled coding sequence: 

LDA 

ENO 

AJP, ZR 
+ LRS 

SCM 

DFAD 



X 

* (if necessary) 

*+3 

48 

=02124000000000000 

=DO0 



IDINT 



PURPOSE: To convert a double precision argument X to fixed point. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($) IDINT, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+3 


00 


($)X 




NOP 






FORTRAN FUNCTION: IDINT (X) 

NORMAL RETURN: The integer result is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN IDINT 

A = (value) Q = (value) 

INTEGER TOO BIG 

This message is written on the standard output unit and a normal exit is taken 
from IDINT when the argument is too large to convert to fixed point. 

STORAGE: 24 locations 

ACCURACY: 35-36 bits 
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INTF 



PURPOSE: To truncate a floating point argument to the nearest floating point integer where magnitude 
is less than or equal to magnitude of the argument. 

COMPASS CALLING SEQUENCE: 



Call by Value (INTF) 




LDA 


X 


L BRTJ 


($)INTF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QINTF) 




L BRTJ 


($)Q8QINTF, ,* 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


00 





FORTRAN FUNCTION: INTF(X) 

NORMAL RETURN: Floating point integer result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 8 locations 

INTF is generally compiled in-line. 

Compiled coding sequence: 

LDA 

ENO 

FAD, UR 



X 

* (if necessary) 

=02044000000000000 
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ITOJ 



PURPOSE: To raise an integer number to an integer power: I**J. 
COMPASS CALLING SEQUENCE: 



Call by Value (ITOJ) 




LDA 


I 


LDQ 


J 


L BRTJ 


($)ITOJ, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name(Q8QITOJ) 




L BRTJ 


($)Q8QITOJ, , * 


L+l SLJ 


*+2 


02 


DICT. 


L+2 00 


($)I 


00 


($)J 



FORTRAN FUNCTION: ITOJ(I,J) 

NORMAL RETURN: The integer result is left in the A register 
ERROR MESSAGE: ERROR DETECTED IN ITOJ 
A = (value) Q = (value) 

Error Code 
Error Codes 

I = 0, J = 0/NEG. 
J GT 47. 
OVERFLOWED INTG. 



I is zero and J is zero or negative 
J is greater than 47 and 1 1 1 > 1 
I I J I is greater than 2 47 -l 



If one of these error conditions occurs, the associated error message is printed 
and a normal exit taken. 



ACCURACY: Exact 

STORAGE: 66 locations including Q2Q07000 

MATHEMATICAL METHOD: 



P*2 +Q+2 1 + R*2 2 + S*2 3 + T*2 4 + U*2 5 



Since J 

(P, Q, R, S, T, U, - .1. or 0) 

J 



Then I 



1 2 "i 4 5 

P*2 Q*2 R*2 S*2 T*2 U*2 
I *j^ *i *i *I *I 



p A Q n c op 

Compute 1,1,1,1 ,1 and multiply the powers of I that are present as indicated by 
the bit pattern of J. 
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ITOX 

PURPOSE: To raise an integer to a floating point power: I**X. 
COMPASS CALLING SEQUENCE: 



Call by value (ITOX) 




LDA 


I 


LDQ 


X 


L BRTJ 


($)ITOX, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by name (Q8QITOX) 




L BRTJ 


($)Q8QITOX, , * 


L+l SLJ 


*+2 


02 


DICT. 


L+2 00 


($)I 


00 


($)X 



FORTRAN FUNCTION: ITOX(I,X) 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: See POWRF for error messages. 

STORAGE: 24 locations including Q2Q07101 

ACCURACY: 34 bits, relative error increases for large arguments 

MATHEMATICAL METHOD: Convert I to floating point and call POWRF. 
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LOGF 



PURPOSE : To compute the natural logarithm of a floating point argument X. 
COMPASS CALLING SEQUENCE: 
Call by Value (LOGF) 



LDA 


X 


L BRTJ 


($)LOGF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QLOGF) 




L BRTJ 


($)Q8QLOGF, , * 


L+l SLj 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 





FORTRAN FUNCTION: LOGF (X) 

NORMAL RETURN: The floating point result is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN LOGF 
A = (value) Q = (value) 

ARG=0/NEG. 

If the argument is zero or negative, this message is written on the standard output 
unit and a normal exit is taken from LOGF. 

STORAGE: 60 locations 

ACCURACY: Exact 75% of the time. Worst known error is 3 in last octal digit. 

TIMING: 154 microseconds 
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LOGF 



MATHEMATICAL METHOD: 1 

X = 2 n f where 2**(-3/8) < f < 2**(5/8) 

LnX = nLn2 + Lnf 

B*Ln2 
Lnf = t (A*Ln2 + + %2 ) - <x * Ln2 

Divide the range into four intervals: 

where 2< 2i_5)/8 < f < 2< 2i_3 ) /8 (for i = 1, 2, 3, 4), 

set = 2 -i / 4 and « = - Ln/8 /Ln2 

Lett= (f -0)/(f +P). 

BLn2 
Then Lnf = t (ALn2 + + ^ ) - a Ln2 

ALn2 = 0. 88924967148818 
BLn2 = -1. 85064969453155 
C = -1.66612572331613 



1 Reference: COOP routine B3 ANL ANLB350 LOGF 
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MAXOF 



PURPOSE: To determine the maximum value of the fixed point arguments L , I„, . . . ,1 . 



COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)MAX0F, , 


L+l 


SLJ 


*+m 




n 


DICT. 


L+2 


00 


<*)Il 




00 


($)I 2 



00 



($)I n 



m = (n+l)/2 + 1, n = number of parameters 

FORTRAN FUNCTION: MAXOF^, i 2 , . . . ) 

NORMAL RETURN: The maximum fixed point argument, converted to floating point is left in the 
A register. 

ERROR MESSAGE: None 

STORAGE: 81 locations for all MIN and MAX functions 



5? 



MAX1F 

PURPOSE: To determine the maximum value of the floating point arguments Xj, X 2 , . . . ,X n . 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)MAX 


:if, , * 


L+l 


SLJ 


*+m 






n 


DICT. 




L+2 


00 


($)X ] _ 






00 


($)X C 





00 ($>x n 

m = (n + 1)/2 + 1, n = number of parameters 
FORTRAN FUNCTION: MAX1F (X p X 2 , . . . ) 

NORMAL RETURN: The largest floating point argument is left in the A register, 
ERROR MESSAGE: None 
STORAGE: 81 locations for all MIN and MAX functions 
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PURPOSE: To determine the minimum value of the fixed point arguments Lp I 2 , . 
the result in floating point. 

COMPASS CALLING SEQUENCE: 



MINOF 



.Ijj, and return 



L 


BRTJ 


($)MIN0F: 


L+l 


SLJ 


*+m 




n 


DICT. 


L+2 


00 


($)Il 




00 


($)I 2 



00 



($)I n 



m = (n+l)/2 + 1, n = number of parameters 
FORTRAN FUNCTION: MINOF (I p I 2> . . . ,1 ) 

NORMAL RETURN: Floating point value of the minimum argument is left in the A register. 
ERROR MESSAGE: None 
STORAGE: 81 locations for MIN and MAX functions 
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MIN1F 



PURPOSE: To determine the minimum value of the floating point arguments X-^ X 2 X n . 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)MIN1F, , * 


L+l 


SLJ 


*+m 




n 


DICT. 


L+2 


00 


($)x x 




00 


($)X 



00 ($)X 

n 

m = (n + l)/2 + 1, n = number of parameters 
FORTRAN FUNCTION: MINIF^, X 2 , . . . ,X n ) 

NORMAL RETURN: The minimum floating point argument is left in the A register. 
ERROR MESSAGE: None 
STORAGE: 81 locations for MIN and MAX functions 
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MODF 



PURPOSE: To compute the value of the floating point argument X.. modulo X„. 
COMPASS CALLING SEQUENCE: 



Call by Value (MODF) 




LDA 


X l 


LDQ 


X 2 


L BRTJ 


($)MODF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QMODF) 




L BRTJ 


Q8QMODF,,* 


L+l SLJ 


*+2 


02 


DICT. 


L+2 00 


($)Xj_ 


00 


($)x 2 



FORTRAN FUNCTION: MODF(X 1 , X 2 ) 

NORMAL RETURN: The value of X 1 modulo X g is left in the A register. If X 1 is zero or the integer 
portion of the quotient Xj/Xg exceeds 2 37 -1, zero is left in the A register. The 
result is negative only if X^ and Xg have unlike signs. 

ERROR MESSAGE: ERROR DETECTED IN MODF 

A = (value) Q = (value) 

Q = ZERO. 

This message is written on the standard output unit if X2=0 and a normal exit 
is taken. 

STORAGE: 38 locations 

ACCURACY: Depends on the size of the arguments and the subtraction, X^ - TX-j/Xgl* Xg. 

MATHEMATICAL METHOD: The result is computed as X x - [x^xj * X 2> where only the integer 

part of the bracketed quantity is used. 



55 



POWRF 



PURPOSE: To perform general exponentiation for A B where A and B are floating point arguments. 
COMPASS CALLING SEQUENCE: 



Call by Value (POWRF) 




LDA 


A 


LDQ 


B 


L BRTJ 


($)POWRF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QPOWRF) 




L BRTJ 


($)Q8QPOWRF, , * 


L+l SLJ 


*+2 


02 


DICT. 


L+2 00 


($)A 


00 


($)B 



FORTRAN FUNCTION: POWRF (A,B) 

NORMAL RETURN: Floating point result is left in the A register 
ERROR MESSAGE: ERROR DETECTED IN POWRF 
A = (value) Q = (value) 

Error Code 
Error Codes 

BASE LT ZERO. 
B*LN(A) GT 709. 
EXP=0/NEG. 



Base is negative 

Maximum is exceeded 

Power is zero or negative and base is zero 



This message is written on the standard output unit and a normal exit is taken 
from POWRF. 

STORAGE: 115 locations 

ACCURACY: 34 bits, as arguments get large relative error increases 

TIMING: 292 microseconds 

METHOD: A = Exponential (B*Ln(A) ) which is computed within the routine. 
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Q1QCPLEX 



PURPOSE: To perform mixed mode conversion for complex mode load and arithmetic operations. 

COMPASS CALLING SEQUENCE: 
Two word operands 



L BRTJ 




($)QlQop3sO 


* 
* > 






L+l DLDA 


($)V,b 




L+2 Return 






One word operands 






L BRTJ 


($)QlQop3sO„* 




L+l ENO 


$V 




LDA 


v,b 




L+2 Return 






Partial word operands 






L BRTJ 


($)Q3Qop340,,* 




L+l 01 


ca,b 




pof 


($)V 




L+2 Return 












V Variable location 


s 


mode of v 




op 


Operation 


b index register 





integer 


00 


load 


ca constant addend 


1 


real 




01 


load complement 


pof parameter offset 


2 


double 




02 


add 




3 


complex 




03 
04 


subtract 
multiply 










Routines included in this pac 


sag 


e 




05 


divide 



Q1Q00300 Load, integer to complex 

Q1Q01300 Load complement, integer to complex 

Q1Q02300 Add, integer to complex 

Q1Q03300 Subtract, integer from complex 

Q1Q04300 Multiply, complex by integer 

Q1Q05300 Divide, complex by integer 

Q1Q00310 Load, real to complex 

Q1Q01310 Load complement, real to complex 
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Q1QCPLEX 

Routines included in this package (continued) 

Q1Q02310 Add, real to complex 

Q1Q03310 Subtract, real from complex 

Q1Q04310 Multiply, complex by real 

Q1Q05310 Divide, complex by real 

Q1Q00320 Load, double to complex 

Q1Q01320 Load complement, double to complex 

Q1Q02320 Add, double to complex 

Q1Q03320 Subtract, double from complex 

Q1Q04320 Multiply^ complex by double 

Q1Q05320 Divide, complex by double 

Q0Q06300 Complement complex AQ register 

Q1Q02330 Add, complex to complex 

Q1Q03330 Subtract, complex from complex 

Q1Q04330 Multiply complex by complex 

Q1Q05330 Divide, complex by complex 

Q3Q00340 Load, logical to complex 

Q3Q01340 Load complement, logical to complex 

Q3Q02340 Add, logical to complex 

Q3Q03340 Subtract, logical from complex 

Q3Q04340 Multiply, complex by logical 

Q3Q05340 Divide, complex by logical 

NORMAL RETURN: Complex result in AQ register 

ERROR MESSAGE: None 

STORAGE: 193 locations 

Q1Q00310, Q1Q01310, and QOQ06300 

are generally compiled in-line. 
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Q1QDUBLE 

PURPOSE: To perform mixed mode conversion for double precision load and arithmetic operations. 
COMPASS CALLING SEQUENCE: 



Full word 






L 


BRTJ 


($)QlQop2sO, ,* 


L+l 


ENO 


$V 




LDA 


V,b 


L+2 


Return 




Partial word 






L 


BRTJ 


($)Q3Qop240, ,* 


L+l 


01 


ca,b 




pof 


($)V 


L+2 


Return 





V Variable location 

b index register 

ca constant addend 

pof parameter offset 



mode of V 



integer 
real 



op 


Operation 


00 


load 


01 


load complement 


02 


add 


03 


subtract 


04 


multiply 


05 


divide 



Routines included in this package 



Q1Q00200 
Q1Q01200 
Q1Q02200 
Q1Q03200 
Q1Q04200 
Q1Q05200 
Q1Q00210 
Q1Q01210 



Load, integer to double 

Load complement, integer to double 

Add, integer to double 

Subtract, integer from double 

Multiply, double by integer 

Divide, double by integer 

Load, real to double 

Load complement, real to double 
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Q1QDUBLE 



Routine included in this package (continued) 



Q1Q02210 
Q1Q03210 
Q1Q04210 
Q1Q05210 
Q3Q00240 
Q3Q01240 
Q3Q02240 
Q3Q03240 
Q3Q04240 
Q3Q05240 
Q0Q06200 



Add, real to double 

Subtract, real from double 

Multiply, double by real 

Divide, double by real 

Load, logical to double 

Load complement, logical to double 

Add, logical to double 

Subtract, logical from double 

Multiply, double by logical 

Divide, double by logical 

Complement double precision AQ register 



NORMAL RETURN: The double precision result is left in the AQ register. 

ERROR MESSAGE: None 

STORAGE: 120 locations 

Q1Q00210, Q1Q01210, and Q0Q06200 

are generally compiled in-line. 
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Q1QREINT 



PURPOSE: To perform conversion from integer to real and from logical to real or integer for 
load and arithmetic operations. 

COMPASS CALLING SEQUENCES: 
Integer to real conversion 



L 


BRTJ 


($)QlQopl00, ,* 


L+l 


ENO 


$V 




LDA 


V,b 


L+2 


Return 




Logical to real or 


integer 




L 


BRTJ 


($)Q3Qopm40, ,* 


L+l 


01 


ca,b 




pof 


($)V 


L+2 


Return 





V variable location 
b index register 
pof parameter offset 
ca constant addend 



m 


mode 



1 


integer 
real 



Routines included 
Q1Q00100 
Q1Q01100 
Q1Q02100 
Q1Q03100 
Q1Q04100 
Q1Q05100 
Q3Q00140 
Q3Q01140 
Q3Q02140 



op 


Operation 


00 


load 


01 


load complement 


02 


add 


03 


subtract 


04 


multiply 


05 


divide 



in this package 

Load, integer to real 

Load complement, integer to real 

Add, integer to real 

Subtract, integer from real 

Multiply, integer by real 

Divide, real by integer 

Load, logical to real 

Load complement, logical to real 

Add, logical to real 
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Q1QREINT 



Routines included in this package (continued) 

Q3Q03140 Subtract, logical from real 

Q3Q 04140 Multiply, real by logical 

Q3Q05140 Divide, real by logical 

Q3Q 00040 Load, logical to integer 

Q3Q01040 Load complement, logical to integer 

Q3Q02040 Add, logical to integer 

Q3Q03040 Subtract, logical from integer 

Q3Q04040 Multiply, integer by logical 

Q3Q05040 Divide, integer by logical 

NORMAL RETURN: The result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 112 locations 

Q1Q00100 and Q1Q01100 

are generally compiled in-line. 
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Q1QSTORE 



PURPOSE: To perform and store mixed mode conversion. 



COMPASS CALLING SEQUENCE: 
Store full word 

L BRTJ 



($)QlQ10msO, ,* 



L+l ENO $V 






STA V,b 






L+2 Return 






Store two full words 






L BRTJ ($)QlQ10msO, ,* 






L+l DSTA ($)V,b 






L+2 Return 






Store partial word 






L BRTJ ($)Q3Q10m40, ,* 






L+l 01 ca,b 






pof ($)V 






L+2 Return 






V Storage location of result 


m,s 


mode 


b index register 





integer 


ca constant addend 


1 


real 


pof parameter offset 


2 


double 




3 


complex 




4 


logical 



Routines included in this package 

Q1Q10100 Store, real to integer 

Q1Q10200 Store, double to integer 

Q1Q10300 Store, complex to integer 

Q1Q10400 Store, logical to integer 

Q1Q10010 Store, integer to real 

Q1Q10210 Store, double to real 

Q1Q10310 Store, complex to real 

Q1Q10410 Store, logical to real 

Q1Q10020 Store, integer to double 
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Q1QSTORE 



Routines included in this package (continued) 

Q1Q10120 Store, real to double 

Q1Q10320 Store, complex to double 

Q1Q10420 Store, logical to double 

Q1Q10030 Store, integer to complex 

Q1Q10130 Store, real to complex 

Q1Q10230 Store, double to complex 

Q1Q10430 Store, logical to complex 

Q3Q10040 Store, integer to logical 

Q3Q10140 Store, real to logical 

Q3Q10240 Store, double to logical 

Q3Q10340 Store, complex to logical 

Q3Q10440 Store, logical to logical 

NORMAL RETURN: The result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 119 locations 

Q1Q10010, Q1Q10310, Q1Q10120, and Q1Q10130 
are generally compiled in-line. 
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Q2Q07000 



PURPOSE: To evaluate I**J. 

COMPASS CODING SEQUENCE: 
L BRTJ 

L+l SLJ 

02 
L+2 LDA 

L+3 LDA 



($)Q2Q07000,,* 

*+3 

DICT. 

($)I 
($)J 



FORTRAN FUNCTION: This routine is called whenever I**J is encountered in a FORTRAN 

statement, and I and J are both type integer. 

NORMAL RETURN: The integer result is left in the A register. 

ERROR MESSAGE: See ITOJ for error messages. 

ACCURACY: Exact 

STORAGE: 66 locations including ITOJ. 

MATHEMATICAL METHOD: See ITOJ. 
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Q2Q07101 



PURPOSE: To evaluate P*X 

COMPASS CODING SEQUENCE: 
L BRTJ 

L+l SLJ 

02 
L+2 LDA 

L+3 LDA 



($)Q2Q07101,,* 

*+3 

DICT. 

($)I 
($)X 



FORTRAN FUNCTION: This routine is called whenever I**X is encountered in a FORTRAN 

statement, and I is type integer and X is type real. 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: See POWRF for error messages. 

STORAGE: 24 locations including ITOX. 

ACCURACY: 34 bits, relative error increases for large amounts. 

MATHEMATICAL METHOD: Convert I to floating point and call POWRF. 



66 



Q2Q07110 



PURPOSE: To evaluate X**I 

COMPASS CODING SEQUENCE: 
L BRTJ 

L+l SLJ 

02 
L+2 LDA 

L+3 LDA 



($)Q2Q07110,,* 

*+3 

DICT. 

($)X 

($)I 



FORTRAN FUNCTION: This routine is called whenever X**I is encountered in a FORTRAN 

statement, and X is type real and I is type integer. 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: See XTOI for error messages. 

STORAGE: 87 locations inch ding XTOI. 

ACCURACY: 36 bits for | 1 1 < 1023, 34 bits for 1 1 1 > 1023 

MATHEMATICAL METHOD: See XTOI. 
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Q2Q07202 

PURPOSE: To compute I**D where I is type integer and D is type double. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q2Q07202, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


LDA 


($)I 


L+3 


DLDA 


($)D 



FORTRAN FUNCTION: This routine is called whenever I**D is encountered in a FORTRAN 

statement. 

NORMAL RETURN: The double precision result is left in the A and Q registers. 

ERROR MESSAGES: Since Q2Q07212 is called, error messages may appear from DLOG and/or DEXP. 

STORAGE: 17 locations 

ACCURACY: 82 bits average 

MATHEMATICAL METHOD: I is converted to floating point and Q2Q07212 (R**D) is called. 
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D POWER 



PURPOSE: To compute D**D, D**R, or R**D where R is type real and D is type double. 

COMPASS CALLING SEQUENCE: 

U = location of base V = location of power 

For D**D 

L BRTJ ($)Q2Q07222, ,* 

L+l SLJ *+3 

02 DICT. 

L+2 DLDA ($)U 

L+3 DLDA ($)V 



For D**R 






L 


BRTJ 


($)Q2Q07221, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


DLDA 


($)U 


L+3 


LDA 


($)V 


For R**D 






L 


BRTJ 


($)Q2Q07212, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


LDA 


($)U 


L+3 


DLDA 


($)V 



FORTRAN FUNCTION: The appropriate Q2Q07. . . routine is called when a U**V is encountered in 

a FORTRAN statement where U and/or V are type double. 

NORMAL RETURN: The double precision result is left in the AQ register. 

STORAGE: 88 (for DEXP) + 124 (for DLOG) + 47 (for DPOWER) = 259 locations 

ACCURACY: 82 bits average 

MATHEMATICAL METHOD: U v is computed as exp(V*LnU), DLOG is called to compute LnU and 

DEXP is called to compute exp(V*LnU). 
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Q2Q07220 



PURPOSE: To compute D**I where D is type double and I is type integer. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q2G 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


DLDA 


($)D 


L+3 


LDA 


($)I 



FORTRAN FUNCTION: This routine is called whenever D**I is encountered in a FORTRAN 

statement. 

NORMAL RETURN: The double precision result is left in the A and Q registers. 

ERROR MESSAGE: ERROR DETECTED IN Q2Q07220 

A = (value) Q = (value) 

D = 0, EXP = 0/NEG. 

This message is written on standard output when the base is zero and the power 
is less than or equal to zero. The most significant part of the base is in the A 
register and the power is in the Q register. 

Since Q2Q07220 calls Q2Q07221 for powers greater than 1023, see Q2Q07221 for 
other error messages. 

STORAGE: 94 locations 

ACCURACY: 83 bits for 1 1 1 < 1023, 82 bits average for |l | > 1023. 

MATHEMATICAL METHOD: The method is the same as that used for XTOI (X 1 ) except the compu- 
tations are in double precision. See METHOD heading of XTOI. 

For powers greater than 1023, I is converted to floating point and 
Q2Q07221 is called. 
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Q2Q07303 



PURPOSE: To compute I**Z where I is type integer and Z is type complex; the real part is in Z; 
the imaginary part in Z + 1. 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q2Q0 


7303, 


* 
> 


L+l 

L+2 
L+3 


SLJ 
02 

LDA 
DLDA 


*+3 
DICT. 

($)I 
($)Z 







FORTRAN FUNCTION: This routine is called whenever I**Z is encountered in a FORTRAN 

statement and I is type integer and Z is type complex. 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See LOGF, POWRF, SINF, and COSF. 

STORAGE: 40 locations including Q2Q07313 

ACCURACY: Depends on POWRF, LOGF, SINF, and COSF and the arguments used. 

MATHEMATICAL METHOD: I is converted to floating point and Q2Q07313 (R** Z) is used. 

(Q2Q07303 and Q2Q07313 are entries in the same routine.) 
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Q2Q07313 



PURPOSE: To compute R**Z when R is type real and Z is type complex; the real part is in Z; 
the imaginary part in Z + 1. 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q2Q07313, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


LDA 


($)R 


L+3 


DLDA 


($)Z 



FORTRAN FUNCTION: This routine is called whenever R**Z is encountered in FORTRAN 

statement and R is type real and Z is type complex. 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See LOGF, POWRF, SINF, and COSF. 

STORAGE: 40 locations including Q2Q07303 

ACCURACY: Depends on POWRF, LOGF, SINF, and COSF and the arguments used. 

MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY X, Y, R are real 

R**Z = R X * cos (Y*Ln(R) )+i*R X *sin(Y*Ln(R) ) 

POWRF, COSF, LOGF, SINF are called to compute R X , 
cos(Y*Ln(R) ), Ln(R), and sin (Y*Ln(R) ) respectively. 
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Q2Q07323 



PURPOSE: To compute D**Z where D is type double and Z is type complex. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($>Q2Q07323, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


DLDA 


($)D 


L+3 


DLDA 


($)Z 



FORTRAN FUNCTION: This routine is called whenever D**Z is encountered in a FORTRAN 

statement and D is type double and Z is type complex. 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See DLOG, EXPF, SINF, and COSF. 

STORAGE: 46 locations 

ACCURACY: Depends on EXPF, SINF, and COSF and the arguments used. 

MATHEMATICAL METHOD: For the complex number Z 

Z = X + iY X, Y are real 

D**Z = D X *cos(Y*Ln(D) ) + i*D X *sin(Y*Ln(D) ) 

where D X = EXP(X*Ln(D) ) 

DLOG, EXPF, SINF, and COSF are called to compute Ln(D), 

EXP(X*Ln(D) ), SIN(Y*Ln(D) ), and COS(Y*Ln(D) ) respectively. 
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Q2Q07330 



PURPOSE: To compute Z**I where Z is complex and I is integer. The real part is in Z; the 
imaginary part in Z + 1 . 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q2Q07330, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


DLDA 


($)Z 


L+3 


LDA 


($)I 



FORTRAN FUNCTION: This routine is called whenever Z**I is encountered in a FORTRAN 

statement and Z is type complex and I is type integer. 

NORMAL RETURN: The real part is left in the A register and the imaginary part in the Q register. 

ERROR MESSAGES: See XTOI, COSF, and SINF. 

STORAGE: 35 locations 

ACCURACY: Depends on XTOI, COSF, SINF, and CABS and the arguments used. 

MATHEMATICAL METHOD: For the complex number Z 

Z = X+IY X and Y are real 

Z**N = |Z| N (COS(N*0) ) +■ i SIN(N*0) ) 
, , 2 2 1/2 

izi = (x + y r^ 

= Tan" (Y/X) such that - tt < 6 < tt 
XTOI, COSF, SINF, CABS, CANG are called to perform 
|Z| N , COS(N*0), SIN(N*0), |Z|, and to get 6, respectively. 



If N is even, compute |z| N = (X 2 + y 2 ) N / 2 
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Q2Q07331.2.3 



Routines Z**R, Z**D, and Z**Z have not been implemented for the complex number Z, the real 
number R, the integer number I, and the double precision number D, because, according to strict 
mathematical definition, these function are either multivalued or indeterminate depending upon 
whether the exponent is rational or irrational. Although computer representation can always be 
expressed as the rational fraction M/N, there are still N legitimate values for the answer. 

To define a unique principal value could lead to frequent discontinuities, and (Z**R)**(l/R) would in 
general be drastically different from (Z** (I/R) )**R. 

Therefore, if one of these routines is called, the following message is written on standard output and 
a normal return is taken: 



ERROR MESSAGE: ERROR DETECTED IN C**S 

A = (value) Q = (value) 

ILLEGAL POWER. 

The values of the A and Q registers are meaningless. 

STORAGE: 11 locations 
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Q2QDLDA 



PURPOSE: To get the arguments for a routine which computes A**B where A and/or B is a two 
word variable (complex or double). 

COMPASS CALLING SEQUENCE: 

The Q2Q type calling sequence is used for all exponentiation, no matter what the types of the 
arguments. 



L 


BRTJ 


($)Q2Q07mst, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


XXX 


($)BASE,b 


L+3 


YYY 


($)POWER,b 



XXX = LDA or DLDA, YYY = LDA or DLDA 



m,s,t 



tv£e 



m = type of result 


integer 


s = type of bas 


e 


1 real 

2 double 


t = type of power 


3 complex 


Then Q2Q07mst may appear: 




Q2Q07mst 


UBJP 


/*\** 


L+l 


XMIT 


(*)*-l, ($)Q8QDICT. 


L+2 


XMIT 


(*)*-2, (*)DICT. 


L+3 


BRTJ 


($)Q2QDLDA, ,* 


L+4 


STA 


GETBASE 




STQ 


GETPOWER 


GETBASE 


OCT 






GETPOWER 



OCT 







NORMAL RETURN: (XXX ($) BASE.b) is left in the A register and (YYY ($) POWER.b) is in the 
Q register. These instructions can then be stored and executed to get the 
actual parameters when needed. 

STORAGE: 7 locations 

USAGE: Q2QLDA is used in all routines processing A**B where A and/or B are type complex 
and/or double. 
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Q2QLOADA 



PURPOSE: To get the arguments A and B to the accumulator for a routine which computes A**B, 
where A and B are single word variables (integer or real). 

COMPASS CALLING SEQUENCE: 

The Q2Q type calling sequence is used for R**R, R**I, I**R, 1**1. I is type integer. 
R is type real. 



L 


BRTJ 


($)Q2Q07mst, ,* 


L+l 


SLJ 


*+3 




02 


DICT. 


L+2 


LDA 


($)BASE,b 


L+3 


LDA 


($) POWER ,b 



where m = type of result 

s = type of base 
t = type of power 

Then Q2Q07mst may appear: 



m,s,t 



1 
2 
3 



tv£e 

integer 
real 
double 
complex 



Q2Q07mst 


UBJP 


A+l 


XMIT 


A+2 


XMT 


A+3 


BRTJ 


A+4 


Renin 



(*)*-l, ($)Q8QDICT. 
(*)*-2, (*)DICT. 
(*)Q2QLOADA, ,* 
Return from Q2QLOADA 

NORMAL RETURN: Upon return from Q2QLOADA, the base will be in the A register and the power 
will be in the Q register. Q2QLOADA executes the actual instructions L + 2, 
L + 3 in the calling sequence to Q2Q07mst to get the base and the power. 

USAGE: Q2QLOADA is used in Q2Q07000 (1**1), Q2Q07111 (R**R), Q2Q07101 (I**R), and 
Q2Q07110 (R**I). 

STORAGE: 10 locations 
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Q7QLODLC 

PURPOSE: To load a logical variable and convert it to floating point. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)Q7QLODLC 

L+l Return 

FORTRAN FUNCTION: This routine is called by a Q3Q. . .routine in a Q1Q. . .package to 

load the logical value specified in a logical mixed mode expression. 

NORMAL RETURN: The floating point result, 0.0 or 1.0, is left in the A register, 

and the Q register is cleared. 

ERROR MESSAGE: none 

STORAGE: 26 locations 
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Q8QDLDA and Q8QDLODA 

PURPOSE: To pick up arguments from successive locations in a Q8Q type of calling sequence. 

COMPASS CALLING SEQUENCE: 

A subroutine with a call by name entry and a two word parameter. (Double, Complex or some 
other 2-word type), with error tracing; may be written as: 



NAME 


UBJP 


/*\** 


L+l 


XMIT 


(*)*-!> ($)Q8QDICT 


L+2 


XMIT 


(*)*-2, (*)DICT. 


L+3 


BRTJ 


($)Q8QDLDA, ,* 


L+4 


Normal return from Q8QDLDA 



where Q8QDICT. is an external symbol. 

The same subroutine with no error tracing may be written as: 

NAME UBJP (*)** 

L+l XMIT (*)*-l, ($)Q8QDCONS 

L+2 BRTJ ($)Q8QDLODA, ,* 

L+3 Normal return for Q8QDLODA 

when Q8QDCONS is an external symbol. 

NAME has the Q8Q type calling sequence. Q8QDLDA uses Q8QDICT. to get the relative 
address of the argument. Q8QDLODA uses Q8QDCONS to get the relative address of the 
argument. 

NORMAL RETURN: The double precision argument is left in the AQ register. 

STORAGE: 14 locations 

USAGE: Q8QDLDA is used to get the double precision argument for DLOG, DSQRT, DSIN/DCOS, 
DEXP. It is also used to get the complex argument for CATAN, CEXP, CSIN/CCOS, and 
CLOG. 

Q8QDLODA is used to get the double precision argument for DCUBRT, and DATAN. It is 
also used to get the complex argument for CANG, CABS, and CSQRT. . 
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Q8QENTRY 



PURPOSE: 1. To initialize storage for the execution of FORTRAN compiled programs by the fol- 
lowing steps: 

Clears Q8QHIST. tables 

Determines and saves available memory 

Writes EXECUTION STARTED AT hhmmss on standard output unit 

2. To terminate input/output operations by calling EXIT. 

Calls IOP. to check and close all input/output tapes 

Returns control to SCOPE 



COMPASS CALLING SEQUENCE: 

Call for initialization after loading a program for execution 



L 


BRTJ 


($)Q8QENTR 


Y,,* 


L+l 


SLJ 


*+l 






00 


DICT. 




Call by EXIT 








L 


BRTJ 


($)EXIT,,* 




L+l 


SLJ 


*+l 






00 


DICT. 





DICT. will De set to contain 
transfer to EXIT 



NORMAL RETURN: 
ERROR MESSAGE: 



Subsequent calls to Q8QENTRY act as a NOP. 
None. 



STORAGE: 51 locations for Tape SCOPE version. 
46 locations for Drum SCOPE version. 
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Q8QERROR 



PURPOSE: To print error diagnostics on the standard output unit. The error diagnostics include: 
Error message 
Contents of A and Q registers 
Name of routine in which error occurred 
A trace of the chain of subprogram calls back to the main program 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q8QERROR, , * 


L+l 


SU 


*+2 




02 


DICT. 


L+2 


00 


($)K 




00 


(DM 



K Location of error key 

M Location of BCD error message 



Alternate entry to preset fatal option 

L BRTJ ($)Q8QERSET, , * 

L+l SLJ *+2 

01 DICT. 

L+2 00 ($)N 

00 



N = 0, all entries to Q8QERROR are considered fatal 

N f 0, all entries to Q8QERROR are treated as specified by (K) 

NORMAL RETURN: If (K) is non-zero, the error message is printed, A, Q registers restored, and 
the next statement is executed. 

FORTRAN FUNCTION: CALL Q8QERSET (N); CALL Q8QERROR (K, M) 

ERROR MESSAGE: (K) = 0, fatal error. The standard printout is followed by: 

EXECUTION DELETED. Control returns to SCOPE. A zero field may 
replace ($)K. 

(K) / 0, non -fatal error. The standard printout is followed by a normal 
return. 

ERROR RETURN: If (K) is zero, job is abandoned. 

STORAGE: 133 locations 
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Q8QIFDIV 



PURPOSE: To interrogate the Interrupt Register to determine if a divide fault, an arithmetic overflow, 
or an exponent overflow has occurred. 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)name, 


L+l 


SLJ 


*H 




00 


DICT. 


L+2 


AJP.MI 


"l 




SLJ 


n 



1 

; name 


condition 


| Q8QIFDJV 

f 

i Q8QIFOVF 
1 Q8QIFEXF 


,. ., \ fault 
divide { , , 
{ check 

arithmetic overflow 

exponent overflow 



FORTRAN FUNCTIONS: IF DIVIDE 



j FAULT / n x 
J CHECK ( 



IF OVERFLOW FAULT n p n c> 
IF EXPONENT FAULT n ± , n 2 

NORMAL RETURN: Control transfers to statement n-^ and the fault condition is cleared if the fault 
has occurred, or to statement n 2 if the fault has not occurred. 

ERROR MESSAGE: None 

STORAGE: 15 locations 
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Q8QIFSSW 



PURPOSE: To interrogate a designated sense switch and branch to one of two given statements, 
depending on the status of the sense switch. 

COMPASS CALLING SEQUENCE: 





ENA 


i if the argument is an integer constant, or 

LDA I if the argument is an integer variable 


L 


BRTJ 


($)Q8QIFSSW, , * 


L+l 


SLJ 


*+l 




00 


DICT. 


L+2 


AJP.MI 


n l 




SLJ 


n 9 



FORTRAN FUNCTION: IF (SENSE SWITCH i) n p n 2 

NORMAL RETURN: If the corresponding sense switch is on, control transfers to statement ny, if the 
corresponding sense switch is off, control transfers to statement n„. 

ERROR MESSAGE: SENSE SWITCH OUT OF RANGE 

If. the integer i is outside the range, 1 ^i^6, 

this message is printed on the standard output unit and 

control transfers to n 2 - 

STORAGE: 20 locations 
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Q8QLOADA and Q8QLODA 

PURPOSE: To pick up an argument. 

COMPASS CALLING SEQUENCE: 

A subroutine with a call by name entry (Q8Qname) and a call by value entry (name), with error 
tracing may be written as: 



Q8Qname 


UBJP 


/*\** 


L+l 


XMTT 


(*)*-l, ($)Q8QDICT 


L+2 


XMIT 


(*)*-2, ($)DICT. 


L+3 


BRTJ 


($)Q8QLOADA 


name 


UBJP 


(*)** 


Lf5 


XMIT 


(*)*-l, ($)Q8QDICT 


L+6 


XMIT 


(*)*-2, ($)DICT. 


L+7 


Normal re 


turn from 080LOADA 



where Q8QDICT. is an external symbol. 
The same subroutine with no error tracing: 



QSQname 


UBJP 


/*\* * 


L+l 


XMIT 


(*)*-l, ($)Q8QLDCON 


L+2 


BRTJ 


($)Q8QLODA 


name 


UBJP 


/*\** 


L+4 


Normal re 


turn from 080LODA 



Q8Qname has the Q8Q type calling sequence. Q8QLOADA uses Q8QDICT. to get the 
relative address of the argument. Q8QLODA uses Q8QLDCON to get the relative 
address of the argument. 

Q8QLOADA returns to *+4 and Q8QLODA returns to *+2. 

NORMAL RETURN: The argument is left in the A register. 

STORAGE: 18 locations 
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Q8QPAUSE 



PURPOSE: To stop operation, displaying an octal constant n on the output comment medium (OCM) . 

1. For PAUSE , a normal return is taken when a character is received from the 
input comment medium (ICM) 

2. For STOPS, a transfer to EXIT in Q8QENTRY is taken when a character is 
received from the ICM. No output occurs on OCM when n is zero. 

COMPASS CALLING SEQUENCE: 

ENA n 

L BRTJ 



($)Q8QPAUSE,,*or ($)Q8QSTOPS, ,* 



L+l 



SLJ 
00 



DICT. 



FORTRAN FUNCTION: PAUSE n 

STOP n 



n may be absent 



NORMAL RETURN: For Q8QPAUSE , operation resumes with the next statement with n in the 

accumulator. For Q8QSTOPS operation resumes, transferring control back 
to EXIT, and the job is terminated. 

ERROR MESSAGE: none 

STORAGE: 32 locations for both Tape SCOPE and Drum SCOPE versions 
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Q8QRESID 



PURPOSE: To initialize actual address and bank parameters in a called subroutine from a parameter 
substitution list.* 

COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)Q8QRESID, , * 


L+l 


SLJ 


Return 




n 


DICT. 


L+2 


r lti 


ca 




Om 


P l 


L+3 


T 2h 


P 2 



n 

m 



77 
77 

number of arguments 
argument ordinal number 



77777,7 
77777,7 



4 




reference 



first reference to argument m 

subsequent reference to same 
argument 



77. 



ca constant addend to address a + Pj 
.7 end of parameter list 



t 


address 





upper address of P= 


1 


upper 18 bits of P. 


2 


upper 24 bits of P; 


3 


lower address of P-, 




upper AUG bank field of P^ 


4 


lower address of P= 


5 


lower 18 bits of P^ 


6 


lower 24 bits of P, 


7 


upper address of P^ + 1, 




lower AUG bank field of P^ 



NORMAL RETURN: The index registers and A,Q registers are restored. 
ERROR MESSAGE: None 
STORAGE: 55 locations 



*Not advisable for handcoded routines. 
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Q8QSENLT 



PURPOSE: 1. To turn off all sense lights, or to turn on a designated sense light. 

2. To interrogate a designated sense light and branch to one of two given statements, 
depending on the status of the sense light. 

COMPASS CALLING SEQUENCE: 
Call by SENSE LIGHT i 

ENA i if i is an integer constant, or LDA I if i is an integer variable 

L BRTJ ($)Q8QSENLT, , * 

L+l SLJ *+l 

00 DICT. 

Call by IF(SENSE LIGHT i) n-p n 2 

ENA i if i is an integer constant, or LDA I if i is an integer variable 

L BRTJ ($)Q8QIFSLT, , * 

L+l SLJ *+l 

00 DICT. 

L+2 AJP.MI nj 

SLJ n„ 



FORTRAN FUNCTIONS: SENSE LIGHT I 

IF (SENSE LIGHT I) n 1; n 2 

NORMAL RETURN: For SENSE LIGHT I: 

If I = 0, all sense lights are turned off; if 1 ^ I < 4, the corresponding sense 
light is turned on. 

For IF (SENSE LIGHT I) n r n 2 : 

If the corresponding sense light is on, it is turned off, and control transfers to 
statement n.; if the corresponding sense light is off, control transfers to 
statement n 2 . 

ERROR MESSAGE: SENSE LIGHT OUT OF RANGE 

If the integer i is outside the range 1 ^ i < 4, this message is printed on the 
standard output unit and a normal return is taken. 

STORAGE: 39 locations 
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Q9QEVAL 



PURPOSE: To compute the offset and operand address of an element in a logical or byte array and 
use this information to position the element in a called subroutine parameter list. 



COMPASS CALLING SEQUENCE: 



L 


BRTJ 


L+l 


n 




f 


L+2 


STA 



name (Return to L+2) 

ca,b 

($)V 

.Z#+i for upper half word or STQ .Z#+i for lower half word 



name 

n 

ca 

b 

f 

($)V 



Q9QEVALL for logical elements 
Q9QEVALB for byte elements 

byte size (number of bits — 01 for logical) 

constant addend 

index register 

bit offset 

base address of array V 



NORMAL RETURN: Upon return, the following information is left in the A and Q registers: 



(A) 

(Q) 

ERROR RETURN: None 
STORAGE: 44 locations 



47 



47 



42, 



,41 39.38 



-t- 
I 
($) I 



(.Z#-+i) 



V + a 



24,23 



I 



(.Z# + i) 



24.23 



18.17 15. 






14 



V+a 
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RANF 

PURPOSE: Repeated use of RANF generates a uniformly distributed sequence of random numbers in 
either fixed or floating point format. If floating point is used, the numbers range from 
to 1; if fixed point is used, the range is to 1 (if interpreted as a fraction) or 1 to 2 47 -1 
(if interpreted as an integer). RANFGET allows extract ion of the current generative 
number, and RANFSET allows it to be replaced. 



COMPASS CALLING SEQUENCE: 
Call by Value (RANF) 

ENA ± 1 
L BRTJ 

L + l SLJ 

00 



($)RANF, ,* 

*+l 

DICT. 



A positive parameter to RANF requests a fixed point response, a negative parameter 
requests floating point. 

Call by Name (RANFGET or RANFSET) 



L 


BRTJ 


($)RANFGET, ,* or (©RANFSET, ,* 


L+l 


SLJ 


♦4-2 




01 


DICT. 


L+2 


00 


($)R 




NOP 






FORTRAN FUNCTIONS: X = RANF(±1) 

CALL RANFGET(R) 
CALL RANFSET(R) 

STORAGE: 33 locations 



ACCURACY: Tests described in COOP routine G5 ANL G550 RANF were performed. 

TIMING: Floatingpoint, approximately 34 microseconds 
Fixed point, approximately 27 microseconds 



MATHEMATICAL METHOD: 



Given the random integer Xj, compute T = 5 15 * X. Then X i+1 = T 
modulo 2 . Since the low order bits of X^ +1 have a very short period, 
the fixed point number returned by RANF consists of bits 47-11 of 
x i+l followed by bits 58-48 of T. The floating point number returned 
is Xj + 1/2 47 . Xj is 11060471625 8 initially. 
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REAL 

PURPOSE: To obtain the real part of a complex number Z. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($)REAL, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)Z 



FORTRAN FUNCTION: REAL(Z) 

NORMAL RETURN: The real part of Z is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 9 locations including AIMAG and CMPLX 

TIMING: 46 microseconds 

REAL is generally compiled in-line. 
Compiled coding sequence: 

LDA Z 
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SIGNF 



PURPOSE: To transfer the sign of argument X 2 to the absolute value of argument Xj. 
COMPASS CALLING SEQUENCE: 

Call by Value (SIGNF or XSIGNF) 





LDA 


X l 




LDQ 


X 2 


L 


BRTJ 


($)SIGNF, , * or ($)XSIGNF, , * 


L+l 


SLJ 


*+l 




00 


DICT. 


Call by Name 






L 


BRTJ 


($)Q8QSIGNF, , * or ($)Q8QXSIGNF, , * 


L+l 


SLJ 


*+2 




02 


DICT. 


L+2 


00 


($)Xi 




00 


($)X 2 



FORTRAN FUNCTION: SIGNFCX^Xg) and XSIGN F(X 1 ,X 2 ) 

NORMAL RETURN: With SIGNF, the sign is transferred between floating point arguments; with 

XSIGNF, the sign is transferred between fixed point arguments. In both cases 
the result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 8 locations 

ACCURACY: Exact 

SIGNF and XSIGNF are generally compiled in-line. 
Compiled coding sequence: 



LDA 


XI 


LDQ 


X2 


QJP.MI 


*+2 


AJP.PL 


*+3 


ROP.XOR 


A,MZ,A 


SLJ 


*+2 


AJP.MI 


*4-l 


ROP.XOR 


A,MZ,A 
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SINF 



PURPOSE: To compute the sine of a floating point argument X in radians. 



COMPASS CALLING SEQUENCE: 




Call By Value (SINF) 




LDA 


X 


L BRTJ 


($)SINF, ,* 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QSINF) 




L BRTJ 


($)Q8QSINF, ,* 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: SINF(X) 

NORMAL RETURN: The value of SINF(X) is in the A register. 

ERROR MESSAGE: ERROR DETECTED IN SIN/COS 
A = (value) Q = (value) 

ARG GT 2**36. 

If the absolute value of the argument exceeds 2 , the message is printed on the 
standard output unit, and the absolute value of the argument is returned in the 
A register. 

STORAGE: 138 locations 

7T 7T 

ACCURACY: Average accuracy in the range ( - — , — ) exceeds 35 bits; the worst known case is 34 bits. 

Li £i 

IT 05 

When |x I > — , the relative error is less than 2 

TIMING: 124 microseconds 

MATHEMATICAL METHOD: 

The method of partitioned polynomials, described in Control Data Technical Report 52 is 
used. Briefly, for |x| less than approximately 50°, a polynomial in x is evaluated; for 
approximately 50° < |x| < 90°, a polynomial in (g - x) is evaluated, in both cases, the 
actual angle determines how many terms in each polynominal are evaluated. 
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SN6L 



PURPOSE: To convert a double precision argument X, to single precision floating point. 



COMPASS CALLING SEQUENCE: 




Call by value (SNGL) 






LDA 


X 


L 


BRTJ 


($)SNGL, ,* 


L+2 


SLJ 


*fl 




00 


DICT. 


Call by name 


(Q8QSNGL) 




L 


BRTJ 


($)Q8QSNGL,,* 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: SNGL (X) 

NORMAL RETURN: Returns with the single precis ion floating point result in the A register. 

ERROR MESSAGE: None 

STORAGE: 24 locations including DBLE 

ACCURACY: 37th bit is rounded into 36th bit. 
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SQRTF 



PURPOSE: To compute the square root of a real argument X. 



COMPASS CALLING SEQUENCE: 




Call by Value (SQRTF) 






LDA 


X 


L 


BRTJ 


($) SQRTF, , * 


L+l 


SLJ 


*+l 




00 


DICT. 


Call by Name (Q8QSQRTF) 




L 


BRTJ 


($)Q8QSQRTF, , * 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






FORTRAN FUNCTION: SQRTF(X) 

NORMAL RETURN: Return with the result in the A register. 

ERROR MESSAGE: ERROR DETECTED IN SQRTF 
A = (value) Q = (value) 

NEG ARG. 

If the argument is negative, this message is printed on the standard output unit and 
a normal exit is taken from SQRTF. 

STORAGE: 31 locations 

ACCURACY: 35 bits average 

TIMING: 107 microseconds 

NOTE: Refer to CO-OP routine B4 ANL ANL B450 SQRTF. 
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SQRTF 



MATHEMATICAL METHOD: 

X = 2 2N * X 1/2 < X< 2 

Vx = 2 N * Vx 

Vx is computed using the Newtonian approximation. 
The initial value is a Chebychev approximation. 1 

B 



Y o 


= A - 


A+X 






A 


= 3.0903156 






B 


= 8.5500505 






Two iterations follow; 








Y i 


= l/2h 


f i-l + 


X 

Yi- 


1 


Then Y 2 « 


-VI 







i = 1, 2 



l W. J. Cody, "Double Precision Square Root for the Control Data Corporation 3600, 
AMD Tech. Memo No. 55, p. 18, Applied Mathematics Division, ANL (unpublished). 
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TANF 



PURPOSE: To compute the tangent of the floating point in radians of X. 

COMPASS CALLING SEQUENCE: 
Call by Value (TANF) 



LDA 


($)X 


L BRTJ 


($)TANF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QTANF) 




L BRTJ 


($)Q8QTANF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


00 





FORTRAN FUNCTION: TANF(X) 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: If the argument is greater than 2 36 : 

ERROR DETECTED IN TANF 

A = (value) Q = (value) 

ARG GT MAX. 
STORAGE: 74 locations 

ACCURACY: TANF (X) was evaluated for 5000 arguments between -5 and +10 radians. The average 
spacing between the arguments was . 003 ; actual spacing was a random number between 
. 002 and . 004. TANF (X) was compared with SINF (X)/COSF (X), which is known to 
have an average accuracy of 1 part in 2 35 , by doing an integer subtraction. 



TIMING: -tt/4 < R < tt/4 

tt/4 < |R| < ?r/2 
tt/2 < |R| 



140 microseconds 
149 microseconds 
187 microseconds 
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TANF 



The differences are as follows: 

Difference -3-2-1 12 3 

Occurrences 202 483 897 1180 869 431 215 
The difference exceeded 3 on 673 samples. Sizable differences occurred only when X 
was approximately equal to k it/2, for k = ±2, ±3, ... 



In this case, when 2 



-n-1 



~r - X < 2 n , an upper bound on the magnitude of the 



relative error is approximately 2]1 . 2 . 



MATHEMATICAL METHOD 



The approximation for computing the tangent of the angle R in radians is given by the 
following equations: 

For -7I-/4 < R < 7r/4 



/ (W) = tan R 



W 



lau t\ - 

S(W Z ) 






W = R * 2/tt 




S(W 2 ) = A + W 2 * 


B 


C 


D + W 2 + E 

F +W2 


A = .63661 


97723 67596 


B = -.07531 


94869 91705 


C = 3.88 


560 


57227 68290 



D =-14.87026 86251 97861 
E =-57.81869 13873 68667 
F = -9.32191 89536 46030 



For I R I > 7r/4, the argument is reduced to its basic equivalent angle by division by w/4. 
Giving R the value of the remainder of this division, the following table shows the values 
of W and /(W). 



Octant 


W 


TANF 


First or fifth 


(2/tt) * R 


/(W) 


Second or sixth 


(2/tt) * (tt/4 -R) 


1//(W) 


Third or seventh 


-(2/tt) * R 


/(W) 


Fourth or eighth 


(2/tt) * (R - tt/4) 


1//(W) 



1 Maehly, Hans J., Approximations for the Control Data 1604 , Control Data Corporation 
publication number 516. 
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TANHF 



PURPOSE: To compute the single precision hyperbolic tangent of a single precision argument X. 

COMPASS CALLING SEQUENCE: 
Call by Value (TANHF) 

LDA X 



L BRTJ 


($)TANHF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QTANHF) 




L BRTJ 


($)Q8QTANHF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: TANHF(X) 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: None 



-34 



STORAGE USED: 71 locations 

ACCURACY: Average 35 bits; worst known relative error 2 

MATHEMATICAL METHOD: 

E 2X = 2 I 2 f = 2 1 *E { * Ln2 

1. TANH(X) = [(E 2X -l.)/(E 2X + l.)J for .3<|X|< 12.48 

q c n Q 

2. TANH(X) = X-— + -^^ - IIX- + 62X_ _ .00886323553X 11 + .00357X 13 

K ' 3 15 315 2835 

for 0< |x| <.3 

3. TANH(0) = 

4. TANH(X) = ±1.0 for |x| > 12.48 
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TIMEF 



PURPOSE: To obtain the current reading of the clock in milliseconds. 
COMPASS CALLING SEQUENCE: 



L 


BRTJ 


($) TIMEF,,* 


L+l 


SLJ 


*+2 




01 


DICT. 


L+2 


00 


($)X 




NOP 






X is a dummy variable 



FORTRAN FUNCTION: TIMEF(X) 

NORMAL RETURN: Floating point result is left in the A register. 

ERROR MESSAGE: None 

STORAGE: 5 locations 
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XDIMF 

PURPOSE: To determine the difference between two fixed point numbers, 1^ and Ig. 
COMPASS CALLING SEQUENCE: 



Call by Value (XDIMF) 




LDA 


h 


LDQ 


h 


L BRTJ 


($)XDIMF, , * 


L-l SLJ 


*+l 


00 


DICT. 



Call by Name (Q8QXDIMF, , *) 



L 


BRTJ 


($)Q8QXDIMF, , * 


L-l 


SLJ 


*+2 




02 


DICT. 


L-2 


00 


($)Il 




00 


($)I 2 



FORTRAN FUNCTION: XDIMF(I y 

NORMAL RETURN: If I, > I„, the fixed point difference is left in the A register. 



If L = !„, zero is left in the A register. 



ERROR MESSAGE: ERROR DETECTED IN XDIMF 
A = (value) Q = (value) 

OVERFLOW ERROR. 



This message is written on the standard output unit and a normal exit is taken 
from XDIMF when (Ij-Ig) I > 2 47 -l. 



STORAGE: 22 locations 

ACCURACY: Exact 

XDIMF is generally compiled in-line. 
Compiled coding sequence: 

LDA XI 

SUB X2 

+ AJP,PL *+l 

ENA 
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XFIXF 



PURPOSE: To convert a floating point number X to the nearest fixed point integer with a magnitude 
less than or equal to the absolute value of X. 

COMPASS CALLING SEQUENCE: 
Call by Value (XFIXF) 

LDA X 



L BRTJ 


($)XFKF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QXFDCF) 




L BRTJ 


($)Q8QXFIXF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: XFIXF(X) 

NORMAL RETURN: Fixed point value is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN XFKF 
A = (value) Q = (value) 

INTEGER TOO BIG. 

This message is printed on the standard output unit and a normal exit is taken 
from XFIXF when the floating point argument is too large to be converted to 
an integer. 

STORAGE: 26 locations 

ACCURACY: Exact 
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XINTF 



PURPOSE: To convert a floating point argument X to the nearest fixed point integer with a 
magnitude less than or equal to the absolute value of X. 

COMPASS CALLING SEQUENCE: 



Call by Value (XINTF) 




LDA 


X 


L BRTJ 


($)XINTF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QXINTF) 




L BRTJ 


($)Q8QXINTF, , * 


L+l SLJ 


*+2 


01 


DICT. 


L+2 00 


($)X 


NOP 






FORTRAN FUNCTION: XTNTF(X) 

NORMAL RETURN: The fixed point result is left in the A register. 

ERROR MESSAGE: ERROR DETECTED IN XFKF 

A = (value) Q = (value) 

INTEGER TOO BIG. 

This message is written on the standard output unit and a normal exit is taken 
from XINTF when X is too large to convert to an integer. 

STORAGE: 26 locations 

ACCURACY: Exact 

NOTE: XINTF is the same routine as XFIXF. 
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XMAXOF 



PURPOSE: To determine the maximum value of the fixed point arguments Ij, I 2 , • • • .Iq. 



IS CALLING SEQUENCE: 




L BRTJ 


($)XMAX0F, , * 


L+l SLJ 


*+m 


n 


DICT. 


L+2 00 


($)Il 


00 


($)I 9 



oo (^ 

m = (n + l)/2 + 1, n = number of parameters 

FORTRAN FUNCTION: XMAXOFfl^ I g I n ) 

NORMAL RETURN: Fixed point value of the largest argument is left in the A register. 

ERROR CODE: None 

STORAGE: 81 locations for MIN and MAX functions 
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XMAX1F 



PURPOSE: To determine the maximum value of the floating point arguments Xj, X2, . . . ,X n . 



COMPASS CALLING SEQUENCE: 




L 


BRTJ 


($)XMAX1F, , * 


L+l 


SLJ 


*+m 




n 


DICT. 


L+2 


00 


($)x x 




00 


($)X 2 



00 



(S)^ 



m = (n + l)/2 + 1, n = number of parameters. 

FORTRAN FUNCTION: XMAXIFpq, X 2 X^ 

NORMAL RETURN: Fixed point value of the largest argument is left in the A register. 

ERROR RETURN: This return is taken if the largest floating point argument results in an absolute 
fixed point value greater than 2 47 -1. 

ERROR CODES: See XFIXF. 

STORAGE: 81 locations for MIN and MAX functions 
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XMINOF 



PURPOSE: To determine the minimum of the fixed point arguments Ij, I 2 . 



■V 



COMPASS CALLING SEQUENCE: 




L 


BRTJ 


($)XMIN0F, , * 


L+l 


SLJ 


*+m 




n 


DICT. 


L+2 


00 


($)Il 




00 


($)I 2 



00 



($)In 



m = (n + l)/2 + 1, n = number of parameters 
FORTRAN FUNCTION: XMINOFflj, Ig, . . . , y 

NORMAL RETURN: Fixed point value of the smallest argument is left in the A register. 
ERROR CODES: None 
STORAGE: 81 locations for MIN and MAX functions 
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XMIN1F 

PURPOSE: To determine the minimum value of the floating point arguments X,, X„ XL. 

COMPASS CALLING SEQUENCE: 

L BRTJ ($)XMIN1F, , * 

L+l SLJ *+m 

n DICT. 

L+2 00 ($)X 1 



OX) ($VX 



o 



00 ($)X n 

m = (n + l)/2 + 1, n = number of parameters 

FORTRAN FUNCTION: XMIN1F(X X . . . ,X ) 

NORMAL RETURN: Fixed point value of the smallest argument is left in the A register. 

ERROR RETURN: This return is taken if the smallest floating point argument results in an 
absolute fixed point value greater than 2 47 -l, 

ERROR CODES: See XFKF. 

STORAGE: 81 locations for MIN and MAX functions 
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XMODF 



PURPOSE: To compute the value of the fixed point argument I modulo J. 
COMPASS CALLING SEQUENCE: 



Call by Value (XMODF) 




LDA 


I 


LDQ 


J 


L BRTJ 


($)XMODF, , * 


L+l SLJ 


*+l 


00 


DICT. 


Call by Name (Q8QXMODF) 




L BRTJ 


($)Q8QMODF, , * 


L+l SLJ 


*+2 


02 


DICT. 


L+2 00 


($>I 


00 


($>J 



FORTRAN FUNCTION: XMODF (I,J) 

NORMAL RETURN: The value of I modulo J is left in the A register. The result is negative only 
if I and J have unlike signs. 

ERROR MESSAGE: ERROR DETECTED IN XMODF 

A = (value) Q = (value) 

DIVISOR IS ZERO. 

This message is printed on the standard output unit and a normal exit is taken 
from XMODF when J is zero. 

STORAGE: 21 locations 

ACCURACY: Exact 
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XTOI 



PURPOSE: To raise a floating point number to an integer power: X**I. 



COMPASS CALLING SEQUENCE: 




Call by 


Value (XTOI) 






LDA 


X 




LDQ 


I 


L 


BRTJ 


($)XTOI, , * 


L+l 


SLJ 


*+l 




00 


DICT. 


Call by Name (Q8QXTOI) 




L 


BRTJ 


($)Q8QXTOI, , * 


L+l 


SLJ 


*+2 




02 


DICT. 


L+2 


00 


($)X 




00 


($)I 



FORTRAN FUNCTION: XTOI(X.I) 

NORMAL RETURN: Return with the result in the A register. 

ERROR MESSAGE: ERROR DETECTED IN XTOI 
A = (value) Q = (value) 

Error Code 
Error Code 

X = 0, I =0 OR NEG. 
EXP OVERFLOW. 



If X is zero and I is zero or negative 
If IX 1 ! ^ 10 308 . 



XTOI calls POWRF for powers greater than 1023. See POWRF for other 
error codes. 



STORAGE: 87 locations including Q2Q07110 
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XTOI 



ACCURACY: 

36 bits for |I| s 1023, 34 bits for |I | > 1023 
MATHEMATICAL METHOD: 

If 1 1 1 > 1023 call POWRF. 

If not, since I = P*2° i-Q+2 1 +R*2 2 +S*2 3 +T*2 4 +U*2 5 +V*2 6 +W*2 ? +Y*2 8 +Z*2 9 
(P,Q,R,S,T ) U,V,W,Y,Z = 1 or 0) 

«. v 1 v P * 2 ° Y ^* 21 V R *2 2 v s *2 3 V T *2 4 v u *2 5 v v *2 6 v w * 2? V Y*2 8 Z*2 £ 
then X=X *X * X * X * X *X *X *X * X *X 

_ . „2 4 v 8 v 16 v 32 v 64 „128 „256 v 512 

Compute X , X , X , X ,X ,X ,X ,X ,X and 

multiply the powers of X that are present as indicated by the bit pattern of I. 
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